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Abstract 

The purely leptonic rare decay Bg fi^ iJi~ is very sensitive to supersymmetric con- 
tributions which are free from the hehcity suppression of its Standard Model diagrams. 
The recent observation of the decay by the LHCb experiment and the first determination 
of its branching fraction motivate a review of their impact on the viable parameter space 
of supersymmetry In this paper we discuss the implications of the present and expected 
future accuracy on BR(Ss for constrained and unconstrained MSSM scenarios, 

in relation to the results from direct SUSY searches and the Higgs data at the LHC. While 
the constraints from BR(S5 can be very important in specific SUSY regions, 

we show that the current result, and even foreseen future improvements in its accuracy, 
will leave a major fraction of the SUSY parameter space, compatible with the results of 
direct searches, unconstrained. We also highlight the complementarity of the Bs ji^ ii~ 
decay with direct SUSY searches. 
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1 Introduction 



The rare decay Bg ii^ iJi~ has been recognised as a probe of new physics beyond the Standard 
Model (SM) and one of the high priority channels for study in the LHC B physics program. 
Because its SM predicted rate is made very small by a helicity suppression, it may reveal the 
contributions of additional diagrams arising in extensions of the SM, which do not suffer from 
the same suppression. In particular, in supersymmetric extensions of the SM (SUSY), its decay 
amplitude receives an enhancement by a factor of order tan^ (3 [l]-[3] , where tan (3 is the ratio 
of vacuum expectation values of the two Higgs fields, and the branching fraction can be larger 
than in the SM by one order of magnitude, or more. The sensitivity of BR(55 1^^ has 
been discussed extensively in the literature in the past years, mostly in constrained versions 
of the Minimal Supersymmetric Standard Model (MSSM) [4-21 and, more recently, in the 



phenomenological MSSM (pMSSM) [22-24 . The recent observation of the decay with the 
determination of its branching fraction by the LHCb experiment to a value very close to the 
SM prediction [25j excludes very large deviations, motivating a review of its implications on 
the viability of SUSY. In this paper, we discuss these implications in the context of constrained 
and unconstrained MSSM models, with R-parity and CP conservation, with an eye also on the 
future experimental progress of this measurement at the LHC. We show that the constraining 
power of BR(55 1^'^ is important in specific regions of the MSSM, but leaves substantial 
room for the SUSY parameters. We quantify this by studying the fraction of the MSSM model 
points, obtained in fiat scans of the parameter space, which are compatible with the present and 
future BR(55 1^'^ constraints in the framework of the Constrained MSSM (CMSSM) and 
the phenomenological MSSM (pMSSM) with 19 free parameters, and account for the results 
on the Higgs mass and direct SUSY searches at ATLAS and CMS. 

The paper is organised as follows. In section 2 we discuss the SM prediction for BR(5s 
11^ the SUSY contributions and the current experimental results. Section 3 discusses the 
experimental prospects at the LHC experiments. The constraints derived are described in 
section 4 for the CMSSM and the more general case of the pMSSM. Conclusions are provided 
in section 5. 



2 Current status 



2.1 SM prediction 

In the SM, the flavour changing neutral current (FCNC) decay Bg ji^ ii~ proceeds via Z 
penguin and box diagrams and is helicity suppressed. The average branching fraction can be 



expressed as 26 -30]: 
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Figure 1: Dominant Bs ji'^iJi' diagrams in the SM, 2HDM and MSSM. 



where /^^ is the decay constant, m^^ is the Bs meson mass and tb^ is its mean hfetime. Cq-^ 
and Cq^ are the Wilson coefficients of the semileptonic scalar and pseudo-scalar operator^ and 
(7io the axial semileptonic Wilson coefficient. The C[ terms correspond to the chirality ffipped 
coefficients. In the SM, only Cio is non- vanishing and it gets its largest contributions from a 
Z penguin top loop (75%) and a W box diagram (24%) (see Fig. [T]). The SM expected value 
is evaluated using mf^{m^) = (4.18 ± 0.03) GeV and mf^ = (173.5 ± 0.6 ± 0.8) GeV (il), 
corresponding to Ciq = —4. 16 ±0.04, from which the following SM prediction for the branching 
fraction is derived 
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BR{Bs /X+/X-) = (3.53 ± 0.38) x 10"^ , (2.2) 
where we used the numerical values of m^, = (5.36677 ± 0.00024) GeV, \VtbVtl\ = 0.0404 ± 



0.0011, tb, = (1.497 ± 0.015) ps pT, 32 and /b, = (234 ± 10) MeV. The value of /b, is 



extracted from the average of the lattice results reported by the ETMC-11 [33], Fermilab- 



MILC-11 [34,35 and HPQCD-12 |36| Collaborations and represents the dominant source of 
systematic uncertainty (8.7%) in the SM prediction. The top mass determination and the choice 
of the renormalisation scheme for its running have an important impact on the evaluation of 
the Bg branching fraction, as discussed in |37|. The effect is illustrated in Fig. [2] where 

we show the SM central value for BR{Bs M^M~) as a function of the top pole mass value. A 
change of ±2 GeV in the top mass corresponds to a ±10~^^ change in the branching fraction 
value. Other sources of uncertainty include the choice of scale for the calculation of the fine- 
structure constant and parametric uncertainties. Adding all these uncertainties in quadrature, 
a total theoretical uncertainty of 11% is estimated. 



2.2 SUSY contributions 

The Bs jJi^ iJi~ decay may receive very large enhancements within specific extensions of the 
SM. In particular, in the MSSM the Higgs-mediated scalar FCNCs do not suffer from the same 

-^Note that Cq-^ ^ = mj) Cs,p- 
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Figure 2: BR{Bs jJi^ ijl ) vs. the top quark pole mass. The black (lower) line corresponds to 
the CP-averaged branching ratio, while the red (upper) line shows the untagged value. 



helicity suppression as the SM diagrams, thus leading to possible drastic enhancements at large 
values of tan /3 yj-[3j. In this case, the Cq-^ , Cq^ coefficients give the dominant contributions. For 
positive values of Cq^ the interference with the term proportional to (Cg^ + Cq^ is destructive. 
The upper bound on BYl{Bs 1^^ is more easily evaded or, conversely, an appropriate 
pseudo-scalar contribution may lead to a suppression of this decay mode to rates below the SM 
expectation. In the MSSM, the largest contribution to Cq-^ and Cq^^ in the large tan^S region. 



reads 1 3, 20 



tan^/5 

Cq, ^ -Cq, ^ (1 + e, tan/5)2 ^ Asm' O^M'^Ml ^^""'^^ ' ^^"^^ 
where x^^ = ji^ ^ with mi the geometric average of the two stop masses, and 

f{x) = --^--^^\nx. (2.4) 

1 — X (1 — xy 

The 65 correction parametrises loop-induced non-holomorphic terms that receive their main 
contributions from higgsino and gluino exchange. Since f{x) > 0, the sign of Cq^ is opposite to 



that of the fiAt term. Here, Eq. (2.3) is given for purely illustrative purposes; in our numerical 
analysis we employ the result of a full calculation, which includes all relevant contributions. It 
must be pointed out that, whereas the MSSM may have a spectacular impact on the Bg 11^ iJi~ 
process, it is equally possible to effectively suppress the SUSY contributions by moving to 
regions of intermediate tan/3 values and/or large masses of the pseudo-scalar Higgs boson 
A. In such cases, the branching fraction does not deviate from its SM prediction, effectively 
preventing this decay from probing parts of the supersymmetric parameter space. 
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2.3 Experimental results 

The Bg ii~ decay has been the target of a dedicated effort at the Tevatron and the 

LHC. To date, the most constraining upper hmit obtained by a single experiment comes from 
LHCb |38|, BR{Bs fi^fi') < 4.5 x 10"^ at 95% C.L., based on 1.0 fb-^ of data at 7 TeV. 



Searches leading to upper limits have been carried out also by CMS 39 and ATLAS |40| , 
while the CDF Collaboration reported an excess of events over the estimated background, 
corresponding to a value BR{Bs = (I-SIq!?) ^ F^l- The combination of the 

LHCb, ATLAS and CMS results led to an upper bound of 4.2 x 10"^ [42] in Summer 2012. 

More recently, the LHCb Collaboration has announced the first evidence for this decay and 
measured its branching fraction p5) to be: 



BR{Bs = (3.2l}:^(stat)igJ(syst)) x 10"^ (2.5) 

This value is in excellent agreement with the SM prediction, leading to speculations on its 
implications on the viability of SUSY. However, it must be noted that the upper limit constraint 
derived from this result is somehow weaker compared to those from the earlier upper limits. 



while it is interesting to investigate the eflFect of the lower limit from (2.5). 

Before discussing these implications, it is important to consider that the theoretical predic- 
tion of the branching fraction does not directly correspond to the quantity measured by the 
LHCb experiment. In fact, the theoretical predictions are CP-averaged quantities in which the 
effect of Bg — Bg oscillations is disregarded. On the contrary, the experimental measurement 
corresponds to an untagged branching fraction which is related to the CP-averaged value by 
the relation 43 , 44 : 

BR{B, ^ /i+/i-)untag = (^^^^) BR(S. ^ IJ,+IJ,-) , (2.6) 

where 

= ^tbAT, = 0.088 ± 0.014 , (2.7) 

and 

_ |Ppcos(2yp)-|^pcos(2yg) 
^Ar - |p|2 + |5|2 ■ (2-8) 

S and P are related to the Wilson coefficients by: 



Cio Ml 1 Cq, - C'q^ 

and 

ips = arg(5) , = arg(P) . (2.11) 
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The resulting untagged branching fraction can be directly compared to the experimental mea- 
surement. The SM expectation for this corrected branching fraction is: 



(3.87 ±0.46) X 10 



1-9 



(2.12) 



In the MSSM, the difference between the CP-averaged and the untagged values of the branching 
fraction depends on the specific SUSY parameters which enter A^r but the shift is typically 
within ±10%. The distribution of the branching fraction values, from our CMSSM scan dis- 
cussed below, is shown in Fig. (S) 



Figure 3: Distribution of BR{Bs M^M~)untag for CMSSM points. The general shape with 
entries at values below and above the SM expected value of 3.87 x 10~^ persists when restricting 
to the points compatible with the results of LHC SUSY searches. 



3 Experimental prospects 

The LHC experiments, in particular LHCb, will keep improving the precision in the deter- 
mination of the Bg ii~ branching fraction. The latest LHCb measurement offers valid 
guidance for estimating the evolution of the measurement accuracy for increasing statistics. By 
symmetrising the statistical uncertainty of the result to L3 and using Gaussian statistics, we 
study the statistical accuracy as a function of the integrated luminosity. At 14 TeV centre of 
mass energy, the production cross section is approximately a factor of two larger compared 
to 7 TeV. The systematic uncertainties are expected to become important once the statistic 
uncertainties drop. These factors are taken into account in the following estimate: 
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for 7 and 8 TeV operations, and 
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(3.14) 
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Figure 4: Expected uncertainty in the branching fraction of Bg ji^ iJi~ vs. the integrated 
luminosity recorded by LHCb (sohd hues). The red (upper) hne refers to an ultimate systematic 
uncertainty of 5% and the green (lower) to an ultimate systematic uncertainty of 1%. The 
dashed lines show the precision of LHC combinations, assuming comparable sensitivity for the 
LHCb and CMS experiments in the same time period. 



for the 14 TeV data, where L is the integrated luminosity, Lq the total integrated luminosity 
taken at 7 and 8 TeV and Ggyst the expected systematic uncertainty. With the improvements 
in the computing power for lattice calculations the uncertainty on /^^ , the dominant source of 
theory uncertainty in the SM prediction, is likely to decrease to ^1% |45|. Fig. [i] shows the 
expected precision in BR(55 1^'^ as a function of the integrated luminosity for LHCb, 
assuming 3.5 fb~^ at 7 and 8 TeV, and two different scenarios for the systematic uncertainties: 
5% and, optimistically, 1%. This shows the importance of improvements in the systematic errors 
over a long period of time. The systematic uncertainty will largely depend on the accuracy 
available for the determination of the fragmentation function ratio fd/ fs- We do not consider 
here improvements to the analysis and the detector performance, which are difficult to quantify 
at present but may lead to an additional reduction of the statistical uncertainties. The upgraded 
LHCb experiment plans to collect 50 fb~^ of data after ten years of running [46], providing an 
ultimate uncertainty of < 2 x 10~^^. In addition, the general purpose experiments can provide 
useful results and the CMS experiment has demonstrated a sensitivity quite close to that of 
LHCb. If this performance can be extrapolated to the future data sets, taking into account the 
larger event pile-up, and the higher energies, the LHC combinations will show improvements 
of > ^/2 on the statistical error compared to the results of LHCb alone. Since the systematic 
uncertainty on fd/ fs is common to all the experiments, it is assumed to be fully correlated in 
this study. 

In summary, we consider two intervals at 95% C.L. for the branching fraction values: 



LI X 10"^ < BR{Bs /X+/X-) < 6.4 x 10"^ 



(3.15) 



corresponding to the current LHCb result of Eq. ( 2. 5] ) and 

3.1 X 10"^ < BR{Bs < 4.6 x 10"^ (3.16) 

which represents a reahstic estimate of the LHC ultimate relative accuracy of ^ 5%, when 
including an estimated improved theory uncertainty of ^ 8% in the determination of the rate 
of this process, if the central value meets the SM prediction. 



4 Constraints in MSSM models 



We study the effect of imposing the constraints of Eqs. (3.15) and (3.16) on the CMSSM and 



pMSSM by performing broad scans of the model parameters and studying the fraction of points 
compatible with those Bs rates. The parameters are varied in flat scans within their 

ranges given below. The SUSY mass spectra are obtained with SOFTSUSY 3.3.4 [47] and the 

We select points where the lightest 



value of BR{Bs with Superlso v3.4 ^,48 



SUSY particle is the Xi neutralino and which are consistent with the LEP and LEP-2 limits 
on SUSY particles. These points are referred to as "accepted" points in the following. Then, 
we test each point for compatibility with the results of the LHC SUSY and Higgs searches. 



4.1 CMSSM 

First, we consider the effect of BR{Bs in the CMSSM parameter space, where we 

perform flat scans varying the CMSSM parameters in the ranges: 

mo,mi/2 G [50,3000] GeV; G [1,60]; G [-10, 10] TeV; sign(/x) > . (4.17) 

Since the effects on BR{Bs are small for negative values of the /x parameter, 

we choose sign(/x) > in the scans, which is also in better agreement with the muon {g — 2) 
constraint. 

Results are given in graphical form in Fig. [5| where we show the BR{Bs values 
as functions of the four CMSSM parameters, comparing the totality of the generated points to 
those consistent with the lightest Higgs boson h mass range, 123 < < 129 GeV [49''50l . The 
BR{Bs admits a lower value of about 1.5 x 10~^, which is still larger than the present 

experimental lower bound derived from the LHCb measurement, so that the experimental lower 
limit does not yet imply the exclusion of portions of the CMSSM parameter space. Branching 
fraction values below ^ 3 x 10~^ can be reached for mi/2 ^ 1 TeV, ^ Aq < 6 TeV and 
tan/5 > 20. However, once the Higgs mass limits are imposed, the vast majority of the allowed 
points have the BR{Bs at values which are equal to, or larger than, the SM prediction, 

with the exception of a few points located in a region at large mo, very small mi/2 ^ 50 — 100 
GeV and Aq ^ 5 TeV. These points are all excluded by e.g. the LEP or Tevatron direct SUSY 
search limits as they lead to too light gluinos and neutralinos. 

As a consequence, in the CMSSM, it is not possible to have BR{Bs smaller than 

the SM prediction and at the same time be in agreement with the SUSY and Higgs search 
results. Therefore, if in the future the central measured value of BR{Bs remains 
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Ao (TeV) tan (3 

Figure 5: Untagged BR{Bs vs. the CMSSM parameters mo (upper left), mi/2 (upper 

right), Aq (lower left), tan^S (lower right). The solid line corresponds to the central value of the 
BR(Bs M^M~) measurement, and the dashed lines to the 2a experimental deviations. The 
green points are those in agreement with the Higgs mass constraint. 



close to the SM prediction, the lower bound is unlikely to have any effect on constraining the 
CMSSM parameter space. 

Figure [6] shows the fraction of CMSSM points compatible with the current LHCb mea- 
surement and the expected ultimate precision in the (mi/2,mo) plane. They are compared 
to the region excluded at 95% C.L. by the ATLAS SUSY searches in channels with missing 
transverse energy (MET) obtained on 5.8 fb~^ of data at 8 TeV [5l] and the expected reach 
of 300 fb~^ at 14 TeV [52j, which shows that the sensitivity through the Bg ii^ ii~ decay 
improves approximately as the reach of direct searches. However, while searches in the jets + 
MET channels are directly sensitive to the mi/2 and mo parameters, the B^ ji^ ii~ decay 
probes a complementary region of the CMSSM parameter space, accessible to direct searches 
only through the H/A tt channel. 

We quantify the fraction of the CMSSM points in agreement with the B>K{Bs 1^'^ 



9 




Figure 6: Fraction of CMSSM points compatible with the current (left) and ultimate (right) 
95% C.L. constraints on BR{Bs M^M~) in the (mi/2,^0) parameter plane. The continuous 
line shows the parameter region excluded by the ATLAS SUSY searches at 8 TeV with 5.8 fb~^ 
of data (from |51|)a nd the dotted line the reach estimated by CMS for searches at 14 TeV with 
300 fb-i (from [52]). 




Figure 7: Fraction of CMSSM points obtained through a 4-parameter flat scan passing the 
LHC SUSY constraints and in agreement with the present BR{Bs M^M~) measurement of 
Eq. (3.15) (continuous line), and with the prospective range of Eq. ( 3.16D (dotted line), as a 
function of tan /3. 



constraint in Fig. [7} As expected the BR{Bs M^M~) provides us with a powerful constraint 
for CMSSM points having large values of tan /3. The fractions of our generated CMSSM points, 
for which we also enforce the requirements to have masses of the sfermions below 3.5 TeV, of 
the gauginos below 3 TeV and of the CP-odd Higgs boson below 2 TeV to make the results 
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Fraction of points 


Current bounds 


Projected bounds 


All CMSSM points 


82.7% 


62.8% 


Accepted CMSSM points 


81.2% 


61.4% 


Points not excluded by LHC searches 


89.2% 


69.0% 



Table 1: Fraction of CMSSM points compatible with the BR{Bs ji^ ji ) constraint. 

directly comparable to those for the pMSSM in the next section, which are compatible with 
the Bs iJi~ rate constraints, are summarised in Table [T| About 11% of the CMSSM 

points not excluded by LHC SUSY searches in our scan are excluded by the current LHCb 
bound. This fraction increases to 31% for the estimated final accuracy on the branching ratio 



of (3.16). We observe that by restricting the analysis to CMSSM points with large values of 
tan/5, i.e. tan/5 > 40, these fractions increase to 21% and 55%, respectively. Instead, imposing 
the anticipated sensitivity of the direct SUSY searches with 300 fb~^ at 14 TeV, the fraction 
of our scan points not excluded by the direct searches and incompatible with the projected 
bounds on Eg decreases from 31% to 23%. 

4.2 pMSSM model 

The pMSSM relaxes the correlations introduced by the mass universality assumptions of the 
CMSSM and allows us to study the inter-relations between the ji^ ii~ yields and the 

MSSM parameters in a general model. Since only a few of these parameters enter in the 
calculation of the Eg fi^fi~ branching fraction, the pMSSM offers also a viable framework to 
study the complementarity of the constraints from this process with those derived from direct 
searches by ATLAS and CMS. 



The analysis performed here adopts the method and tools described in [22,53 . We perform 
fiat scans of the 19 pMSSM parameters in the ranges: 

Mi,M2 G [-2500,2500] GeV; M3 G [50,2500] GeV; tan/5 G [1,60] 

Ma G [50, 2000] GeV; A, A, G [-10, 10] TeV; /x G [-3, 3] TeV (4.18) 

^Ilr ^ [50,2500] GeV; rriq^ ^ G [50,3500] GeV . 

The dependence of the BR{Bs values calculated at each pMSSM point with the 

most relevant pMSSM parameters is given in Fig. |8] for all the valid points and those having 
123 < Mh < 129 GeV. Contrary to the case of the CMSSM, here even after imposing the 
Higgs mass constraints a sizeable number of points with a value of BR(Ss below the 

SM prediction (down to 0.5 x 10~^) is obtained. These low values are reached for tan/5 > 10 
and > 300 GeV. This observation is important for the prospect of improving the lower 
experimental bound on the decay ratej^ 

^Note that the lower reachable value we obtain is smaller than the one obtained in the recent study of 
Ref. This is because we use the full MSSM expressions with no assumption on Cio, and use non universal 
masses for the SUSY particles. 
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Figure 8: Untagged BR{Bs — )• fJ-'^lJ^~) vs. the parameters n (upper left), M3 (upper right), 
At (middle left), tan/3 (middle right). Ma (lower left) and (lower right). The solid line 
corresponds to the central value of the BR{Bs — )• yU'''yU~) measurement, and the dashed lines to 
the 2a experimental deviations. The green points are those in agreement with the Higgs mass 
constraint. 
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^10 

Figure 9: Variation of the untagged BR{Bs in the plane (Cio,CgJ. The dotted 

vertical lines delimit the range of Ciq in the CMSSM, and dashed lines the range in the pMSSM. 




Figure 10: Constraints from BR{Bs in the (M^,tan/5) and {Ma^tti^J parameter 

planes. The black points corresponds to all the valid pMSSM points and those in grey to the 
points for which 123 < < 129 GeV. The dark green points in addition are in agreement 
with the latest BR(5s range given in Eq. ( |3.15| ), while the light green points are in 

agreement with the prospective LHCb BR{Bs range given in Eq. (3.16). The red line 

indicates the region excluded at 95% C.L. by the CMS A/H r^r~ searches (from [54j). 



The BR(Ss 1^'^ dependence on the Cio and Cq-^ = —Cq^ Wilson coefficients in the 
minimal flavour violation (MFV) framework [55, 56 is shown in Fig. |9| It is instructive to 
observe that the values of BR(Ss 1^'^ can decrease down to for Cio = Cq-^ = 0. 
However, in the pMSSM, the variation of Cio is limited to the interval [-5.0,-2.6], even when 
applying constraints from B observables, so that the lowest value which can be 
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Figure 11: Fraction of pMSSM points passing the LHC SUSY and Higgs mass constraints and 
in agreement with the latest BR{Bs measurement of Eq. (3.15) (continuous hue), 



and with the prospective range of Eq. (3.16) (dotted hue), as a function of tan/5. 



reached is around 0.5 x 10 ^0 

The impact of the present and future determinations of BR{Bs on the parameters 



most sensitive to its rate: (M^, tan /5) and (M^, J is shown in Fig. 10, where we give aU the 



vahd pMSSM points from our scan, those with 123 < < 129 GeV and, highhghted in green. 



those in agreement with the present BK{Bs /x^// ) range (3.15) and the ultimate constraint 
(3.16) at 95% C.L. As already discussed in Ref. [22], the constraints from BR{Bs 



affect the same pMSSM region, at large values of tan/3 and small values of M^, also probed 
by the dark matter direct detection constraints and, more importantly, the H/A T^r~ 
direct Higgs searches at the LHC |54|[58]. The search for the H/A t^t~ decay has already 
excluded a significant portion of the parameter space where large effects on BR{Bs 
are expected. We also note that the stop sector is further constrained by direct searches in 
6-jets + MET channels, which disfavour small values of m^^. The figure shows that it is difficult 
for Ma and m^^ to be simultaneously light. 

In more quantitative terms, we compute the fractions of all the accepted pMSSM points, 
of those not excluded by the jets + MET and H/A t^t~ searches by ATLAS ^58j and 
CMS [54], and those also compatible at 90% C.L. with th e ATL AS a nd CM S Higgs data using 
the analysis of Ref. |59|, which are compatible with the (3.15) and (3.16) constraints at 95% 
C.L. on Bs jji^ iJi~ . Results are summarised in Fig. 11 and Table [2} The current LHCb result 
rules out just below 3% of the pMSSM points compatible with the LHC direct SUSY searches 
and the Higgs results. The projected bound, assuming the central value coincides with the SM 
expectation, will increase the reach by an order of magnitude to 30% of the points and severely 
constrain solutions with very large values of tan/3. By then, the direct searches for SUSY in 
channels with jets + MET and H/A t^t~ will also have extended their sensitivity to a much 
larger part of the pMSSM parameter space. Extrapolating the current bounds to 300 fb~^, the 
fraction of our scan points not excluded by the direct searches but excluded by the projected 



^In general non-SUSY MFV scenarios, Cio can admit larger ranges, leading to constraints also coming from 
the lower bound of BR(55 as shown in ^57j. 
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Fraction of points 


Current bounds 


Projected bounds 


All pMSSM points 


95.3% 


67.8% 


Accepted pMSSM points 


97.7% 


78.1% 


Points not excluded by LHC searches 


95.1% 


63.3% 


Points compatible at 90% C.L. with Higgs results 


97.2% 


70.0% 



Table 2: Fraction of pMSSM points, obtained through a 19-parameter flat scan, compatible 
with the BR(55 constraint. 

bounds on Bg ji^ ii~ will decrease from 30% to c::^20%. 

If SUSY is indeed realised in nature and a signal from the direct searches at ATLAS and 
CMS is observed by then, it would be interesting to perform a quantitative test of consistency 
between the mass of the SUSY states being observed and the branching fraction of this, until not 
long ago, elusive decay. In particular, if the pseudoscalar Higgs and the scalar top masses are 
determined by ATLAS and CMS, the precise value of BR(5s 1^'^ can be used to severely 
constrain the combination of (/i^t, tan/3) in the MSSM. Moreover, by constructing alternative 
observables such as double ratios of leptonic decays formed from the decays 
Bu D ^ jiv and Ds jiv/rv^ it could be possible to enhance the sensitivity of the 

individual decays, through the cancellation of hadronic uncertainties, and stronger constraints 
can be obtained [13^,16] . 

5 Conclusions 

The observation of the rare decay Bs ji^ iJi~ and the flrst determination of its branching 
fraction by the LHCb experiment represent a major milestone of the probe of physics beyond 
the SM through rare decays of h hadrons. The excellent agreement of the measured value 
with the SM prediction has raised the question of its implications on the viability of SUSY. 
In this paper, we have reviewed the predictions for the branching fraction of this decay in the 
SM and the MSSM and discussed the impact of the new LHCb result and the expected flnal 
LHC accuracy on the SUSY parameter space in two models: the CMSSM and the pMSSM. 
We observe that, despite the signiflcant differences between the two models, the sensitivity of 
the Bs ji^ ii~ rate is signiflcant in speciflc regions of the parameter space of these models, 
mostly at large values of tan/3, regions which are also probed by direct SUSY particle searches 
at ATLAS and CMS. As a result, the constraint derived from the current LHCb result removes 
^10% of the scan points in the CMSSM and a few % in the pMSSM, which are not already 
excluded when the bounds from direct SUSY searches and the Higgs data are applied. This is a 
consequence of the suppression of the SUSY contributions for intermediate tan/3 values and/or 
large masses of the pseudo-scalar Higgs boson ^4, where the branching fraction in the MSSM 
does not deviate from its SM prediction. The situation in other constrained MSSM scenarios, 
such as AMSB and GMSB, is similar to that in the CMSSM, with high sensitivity at large 
tan /3 flG^ . The improved accuracy on the branching fraction measurement expected from the 
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14 TeV runs together with the expected improvements in the theory uncertainties, will boost the 
sensitivity, in particular for the region tan > 50 which could be almost completely constrained, 
and underline the complementarity of direct and indirect searches for supersymmetry through 
the possibility of consistency checks, if the heavy Higgs bosons can be observed in the direct 
searches conducted by ATLAS and CMS. 
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